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Abstract 

We study the error induced by the time discretization of a decoupled forward- 
backward stochastic differential equations (X, Y, Z) . The forward component X is 
the solution of a Brownian stochastic differential equation and is approximated by 
a Euler scheme X N with N time steps. The backward component is approximated 
by a backward scheme. Firstly, we prove that the errors (Y N -Y,Z N - Z) measured 
in the strong L p -sense (p > 1) are of order N~ 1 / 2 (this generalizes the results by 
Zhang [20]). Secondly, an error expansion is derived: surprisingly, the first term is 
proportional to X N — X while residual terms are of order iV -1 . 

Key words: backward stochastic differential equation, discretization scheme, 
Malliavin calculus, semi-linear parabolic PDE. 
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1 Introduction 



Let (O, J 7 , P) be a given probability space on which is denned a g-dimensional 
standard Brownian motion W, whose natural nitration, augmented with P- 
null sets, is denoted by (J r t)o<t<T (T is a fixed terminal time). We consider 
the solution (X, Y, Z) to a decoupled forward-backward stochastic differential 
equation (FBSDE in short). Namely, X is the Revalued process solution of 



X t = x+ [ b(s,X s )ds+ f* cr(s,X s )dW s 
Jo Jo 
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and Y (resp. Z) is a real-valued adapted (resp. predictable I^-valued) process 
solution of 



-dY t = f(t,X t ,Y t ,Z t )dt-Z t dW t , y r = $(x r ). (2) 

We assume standard Lipschitz properties on the coefficients, which ensure 
existence and uniqueness in appropriate L 2 -spaces (see Pardoux and Peng 



[18j . or Ma and Yong [14j for numerous references). During the last decade, 
more and more attention has been paid to these equations, because of their 
natural applications in Mathematical Finance or in the probabilistic resolution 
of semi-linear partial differential equations (PDE in short): see El Karoui et al. 
[11 or Pardoux 17 . 



Our aim is to study the most usual time approximation of (X, Y, Z). For X, 
we use the Euler scheme X N with N discretization times (t k = kh) <k<N 
(h = £ is the time step). For convenience, set AW k = W tk+1 — W tk (AW^ 
component- wise) . X N is defined by Xq = x and 

t e [h,t k+1 ], X t N = x t N k + b(t k ,X t N k )(t - t k ) + *(t h ,X»)(W t - W tk ). (3) 

The backward SDE (2) is approximated by (Y N , Z N ) defined in a backward 
manner by Y t N = $(X^) and 

Y t N k =^ fc (0 + hE t J(t k , X», Y£ +i ,Z»), (4) 
hZ» =E tk (Y» + AWZ), (5) 

where K tk is the conditional expectation w.r.t. Tt k and * is the transpose 
operator. Additional tools are needed to derive a fully implementable scheme, 
in particular for the computations of conditional expectations. We refer to 
Bouchard and Touzi for Malliavin calculus techniques, or to Gobet et al. 
0, Lemor et al. 0] for empirical regression methods. In this work, we leave 
these further questions and we only address the error analysis between (Y, Z) 
and (Y N , Z N ). 

On the one hand, Zhang 2oj proves (in a slightly different form) that the error 



maxfc<7v || Y t N — Y tk \\l 2 < CN^ 1 ^ 2 . This is done under rather minimal Lipschitz 
assumptions on b, a, f, $. On the other hand, when / does not depend on 
z and the coefficients are smooth, one knows that \Y^ — Y Q \ < CN~ X (see 
Chevance jjj). We aim at filling the gap regarding these two different rates of 
convergence. In the following, we prove that 

• the Chevance's results extend to the case of / depending also on z. 

• the rate N~ l holds true also for the difference \Zq — Z \. 

• more generally, for the other discretization times t k , we expand the error as 

Y£ - Y tk - a k ■ (X t N k - X tk )\ < CN~ l V \X» - X t f 
(for an explicit and bounded random vector ot k ). 
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• an analogous expansion is available for Z. 

Since \Xf — X tk \ 2 has the same order in L p than N~ l , the error on Y is mainly 
due to the error Xj* — X tk . Thus, Zhang's results are a consequence of this 
expansion, and Chevance's ones as well since Xq = Xq. The gap is filled. 
In addition, we learn from this expansion that if one could perfectly simulate 
X (as for Brownian motion with constant drift, geometric Brownian motion 
or Ornstein-Uhlenbeck process), the error on the BSDE would be of order 
N^ 1 and not N~ l l 2 as stated by Zhang's results. Also, if one could use a dis- 
cretization scheme for X of order 1 for the strong error (for instance Milshtein 
scheme whenever possible), the error on the BSDE would be of order N~ x 
(we would need to extend our analysis to other discretization schemes, this is 
straightforward for the Milshtein scheme). 

The paper is organized as follows. In Section 2, we define the assumptions on 
the coefficients, recall the connection between BSDEs and semi-linear PDEs 
(which is important for our analysis). Finally, we state our main results. Firstly 
in Theorem 6, we extend Zhang's results to L p norm. Secondly in Theorem 7, 
we expand the error on Y. Lastly in Theorem 8, we deal with the error on Z. 
Naturally, stronger and stronger assumptions are required for theses theorems. 
Proofs of the three results are postponed to Sections 3, 4 and 5: we combine 
BSDE techniques, martingale estimates and Malliavin calculus. 

Notation. 

• Differentiation. If g : R d i— ► R q is a differentiable function, its gradient 
V x g(x) = (d Xl g(x), d Xd g(x)) takes values in R q ® R d . At many places, 
V x g(x) will simply be denoted g'(x). If g : R d i— > R is a twice differentiable 
function, its Hessian H x (g) takes values in R d (g) R d : (H x (g))i t j = d 2 x . x ^g. If 

g : R d x R q h+ R, £ y (x,y) takes values in R d ® R q : {£ y )ij = ^Q^, for 

1 < i < d,l < j < q. 

• Function spaces. For an integer k > 1, we denote by (j^ 2 ' k ' k ' k the set of 
continuously differentiable functions <ft '■ (t, x, y, z) G [0, T] x R d x R x R q i— > 
<f)(t,x,y,z) such that the partial derivatives d^d^d^d^^t^x^y^z) exist for 
21q + li + l 2 + h < k and are uniformly bounded. The analogous set of 
functions that not depend on y and z is denoted by C^ 2,k . This set is 
denoted by (j( k+a )l 2 > k+a ( a e]0, 1[) if in addition the highest derivatives 
are Holder continuous with index a w.r.t. x and a/2 w.r.t. t (for a precise 
definition, see Ladyzenskaja et al. |l2|). 

• Norm. For a ci-dimensional vector U, we set \U\ 2 = J2i=iU 2 . For adx 
g-dimensional matrix A, A t denotes its z-th column, and A 1 its z-th row. 
Moreover, \A\ 2 = £j? = iA ? j- 

• Constants. Let C denote a generic constant which may depend on the 
coefficients b, a, f, $ and on the dimensions d and q. We will keep the same 
notation K{T) for all finite, nonnegative, and nondecreasing functions w.r.t. 
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T: they do not depend on x and h. The generic notation K(T, x) stands for 
any function bounded by K(T)(1 + \x\ q ), for some q > 0. 

• 0(U) and Ok(h). A random vector R is such that R = 0(U) for a nonneg- 
ative random variable U if \R\ < K(T,x)U (in particular, i? = 0(/i) means 

< K(T, x)h). The notation R = Ok(h p ) means \R\ < h p , where is 
Tt k -measurable, sup Ar E(sup fc (A^l 9 ) < K(T,x), for q > 1. 

• E tfe and Var tfe . K tk is the conditional expectation w.r.t. Tt k and Var( fc (X) = 
E tk (X*)-(E tk (X)f. 

• Malliavin calculus. We use the notations of Nualart [l6( for weak spaces 
B k ' p . 

• Discretization Let s G [t&,£fc+i[- We define r](s) = 



2 Main results 



2. 1 Hypotheses 



The coefficients b : [0, T]xR d ^ R d , a : [0, T] x R d -> M dx ", / : [0, T]xK J x 
M x K ? ^ R and $ : R d —>■ R satisfy one of the following set of assumptions. 

Hypothesis 1 The functions b, a, f and $ are bounded in x, are uniformly 
Lipschitz continuous w.r.t. (x,y,z) and Holder continuous of parameter | 
w.r.t. t. In addition, $ is of class C% +a for some a g]0, 1[ and the matrix- 
valued function a = aa* is uniformly elliptic. 

2 3 

Hypothesis 2 Assume Hypothesis 1 and that the functions b, a are in C 6 2 ' , 
f is in Cff' 3 ' 3,3 , $ is in C^ +a for some a e]0, 1[. 

Hypothesis 3 Assume Hypothesis 1 and that the functions b, a are in C b ' , 

f is in C^' 4,4 ' 41 , $ is in C A+a for some a g]0, 1[. 

We do not assert that these smoothness and boundedness conditions are the 
weakest ones for our error analysis, but they are sufficient. Investigations re- 
garding minimal assumptions would be certainly interesting but it is beyond 
the scope of the paper. 
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2.2 Connection between Markovian BSDE's and semi-linear parabolic PDE's 



We recall classical results connecting (Y, Z) and the solution and its gradient 
of the following semi-linear PDE on [0, T] x R d : 

(d t + C(t, x ))u(t, x) + / (t, x, u(t, x), V x u(t, x)a(t, x)) = 0, (6) 
u(T,x) = $(z), 

where C{t, x ) is the second order differential operator 

(see for instance Ma and Zhang 3| or Pardoux 
Proposition 4 Under Hypothesis 1, one has 

Vt G [0,T], y t = u(t,X t ), Z t = V x u(t,X t )a(t,X t ), (7) 
where u is the unique classic solution Cl' 2 of the PDE (6). 

- 3 2 4 

In addition under Hypothesis 2, u G C 6 2 ' , and under Hypothesis 3, u G C b ' . 



The first result of this Proposition corresponds to Theorem 2.1 of Delarue and 
Menozzi The two last regularity results can be proved in the same way. 
In fact for this, we would only need 6, a to be in (j^ +a l 2 > 2+a ^ fae additional 
smoothness is used later for Malliavin calculus computations. 



2.3 Main results 



We now turn to the statement of our results. Remind the following well-known 
upper bound on the Euler Scheme, which is useful in the sequel. 

Proposition 5 Let o and b be Lipschitz continuous. Then 
Wp> l,[E(sup|Xf -X t \ p )}» <K(T,x)^=. 



In fact, for all p > 1 one has 

[E^( sup \X t N -Xtffp <K(T,X ti )^= + \X»-X ti \. (8) 

U<t<T \JN 

Our first result is an extension of the L 2 estimates in Zhang [2(j to L q estimates 
(see also Gobet et al. 0). 
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Theorem 6 Let us assume Hypothesis 1. Let q > 0. We define the error 

N-l . t i 

e q (N) = [maxE|y tfc -Y»\* + EC£, / |Z£ - ^| 2 di)*]% 

/c=0 fc 

where Y N and Z N are defined by (4) and (5). Then \e q {N)\ < K(T,x)^=. 

By slightly strengthening the smoothness assumptions on b, a, f and $, we are 
able to expand the error on Y . 

Theorem 7 Let us assume Hypothesis 2. Then, the following expansion holds 
Y t N k - Y tk =V x u(tk,X tk )(X t N k - X tk ) + O k (h + 0(\X t N k - X t f ). 



In view of Proposition 5, \XF — X tk \ 2 and A^ -1 have the same order (in L p ). 
Hence it turns out that V x u(tk, X tk )(X^ — X tk ) is the first order term in 
the error Y t ^ — Y tk . Obviously, this estimate implies that of Theorem 6. As 
mentioned in the introduction, the evaluation of Y by Y^ has still an accuracy 
of order A^ -1 since initial conditions for X N and X coincide. Note that if there 
is no discretization error for the process X, Y t ^ — Y tk = O(jj), a fact which 
is not clear from equations (4) and (5). A nice situation corresponds to a 
independent of x (this is a very specific situation where Euler and Milshtein 
schemes are equal): in that case \\X% - X tk \\ Lp = OiN- 1 ) and one gets the 
order of accuracy A^ 1 for Y. 

For Z which plays the role of a gradient relatively to Y, we get an analogous 
result about the error, up to increasing by 1 the degree of smoothness of the 
coefficients. 



Theorem 8 Let us assume Hypothesis 3. Then, the following expansion holds 
Z» - Z tk = (V X [V X « a]*(t k , X tk )(X t N k - X tk ))* + O fc (-L) + (\X t N k - X tk | 2 ). 

Remark 9 The above results are sufficient to derive the weak convergence of 
the renormalized error process [\^N(Y t N — ^t)]o<t<T an d [^/^{Z^ — Z t )] <f<T; 
except that one has to define Y N and Z N between discretization times. For 
t G [tk,tk+i[, analogously to (4) and (5) we define 
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Theorems 7 and 8 can be extended to all t £ [0,T]. We have 

Y t N - Y t =V x u(t, X t )(X t N - X t ) + O t (i) + 0(|Xf - X| 2 ), 

Zf - Z t =(v,[V,ti <x]*(t, X t )(Xf - X,))* + Ot(i) + 0(|Xf - X| 2 ). 

Theorem 3.5 of Kurtz and Protter fUlJ allows us to establish the weak conver- 
gence of the processes \/N(Y N — Y), and y/N(Z N — Z). Indeed, the process 
[a/X(X/ v — X t )] <t<T weakly converges to the solution of 

Ut=fl [ t V x a l (s,X s )U s dW l s + fVsbfaXjUsds 
i=1 Jo Jo 

i,j=i jQ k=l 

where (V l ^)x<i,j<q ore independent standard Brownian motions and indepen- 
dent ofW. Furthermore, the convergence is stable (see Jacod and Protter 
Hence, [VN(X t N -X t ), y/N (Y t N —Y t ), \/N (Z? - Z t ) , X t ] Q < t < T weakly converges 
to [{U t ,V x u{t,X t )U t ,{[V x [V x u <j}*{t,X t )}U t )*,X t } <t<T. 



3 Proof of theorem 6 

Extra notations for all the proofs. For any process U (except the Brownian 
increments AWk), we define AUk = — U tk . Let 9 S denote (s, X s , Y s , Z s ) 
and /£ denote Y^.Z^). 

Z tk is defined as hZ tk := E tk Z s ds and we put AZ fc = — Z tk . 

If q = 2, the result has already been proved in Gobet et al. p, under Lipschitz 
conditions on b, a, f, $. Thanks to the inequality E\U\ q < (E\U\ 2p )^ for 2p > 
q, we only need to prove the theorem for q = 2p, where p EN*. 
First, we give some estimates which can be easily established. We have, under 
Hypothesis 1, Vs G [t k ,t k+1 ], 

E tk (\X s - X tk \* + \Y S - Y t f* + \Z S -Z tk \ 2 n < ChP. (9) 

In the following computations, these estimates are repeatedly used. 

3. 1 Proof of max <k<N E| Y th - Y t N k \ 2 ? = O (hP) . 

We prove the following result, which is a bit more general. 

Proposition 10 max^^vE^ - Y t f \ 2 ? = O^hP) + |AX| 2 ^. 
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By taking % = 0, we get maxo< fc <ivE|Y tt - Y t f \ 2p = 0(h p ). 
Assume that we have 

\AY k \ 2 < (1 + Ch)E tk \AY k+1 \ 2 + Ch\AX k \ 2 + Ch 2 . (10) 

Then, using the inequality (a + b) p < a p (l + e(2 p ' 1 - 1)) + b p {l + ^r 1 ) for 
< e < 1, we deduce 

\AY k \ 2p < (1 + Ch) p+1 E tk \AY k+1 \ 2p + C p h p (\AX k \ 2 + Ch) p (l + 

Take the conditional expectation w.r.t. T u to get E ti | AY k \ 2v < {l+Ch)E u \ AY k+1 \ 2p + 
h(h p +E t . | AX k \ 2p ). Using (8) for | AX k \ and Gronwall's lemma yields m&Xi< k < N E t . \Y tk 
Y t ^\ 2p = Oi(h p ) + \AXi\ 2p . □ 

Now we prove the inequality (10). From (2) and (4) we obtain 

AY k = E tk (AY k+1 )+E tk [ tk+ \f t N k - f(9 s ))ds. (11) 

By applying Young's inequality, that is (a + b) 2 < (1 + 7/1) a 2 + (1 + ^)b 2 , 
where 7 will be fixed later, and using the Lipschitz property of /, we get 

|AY fc | 2 < (1 + 7 / i )(E tfc (AY fc+1 )) 2 + C(h + -)[h 2 + E tk / \X S - X t N k \ 2 ds] 

7 ^fc 

+ C(h+k[E tk f k+1 \Y s -Y t N k J 2 ds + E tk \ th+1 \Z s -Zl\ 2 ds\. (12) 

7 Jt k J tk 

Let us introduce Z tk (see extra notations at the beginning of Section 3): 

|Z. - ^| 2 rf S = E tk [ tk+1 \Z s -Z tk \ 2 ds + hE tk \Z tk - Zf k \ 2 . (13) 
Jt k Jt k 

Thanks to the Cauchy Schwarz inequality we have 

|E tfc (AY fc+1 A^)| 2 < h{E tk (\AY k+1 \ 2 ) - |E tfc (AY fc+1 )| 2 }. 

Hence, as hZ tk = E tfc ({Y tfe+1 + J^ +1 f(9 s )ds}AW* k ), with a bounded /, it 
follows that 

h 2 \Z tk -Zl\ 2 <dh (E tfc (|AY fc+1 | 2 ) - |E tfc (AY fc+1 )| 2 ) + Ch\ (14) 

By plugging (13) and (14) into (12), we get 

|AY fe | 2 < (l + 7 /i)(E tfc (AY fe+1 )) 2 

+ C(h + -)[h 2 + E tk f k+1 \X s -X t N k \ 2 ds + E tk f k+1 \Y s -Y t N k J 2 ds} 

7 Jt k J tk 

+ C(h + -)[E tk f k+1 \Z s -Z tk \ 2 ds + E tk {\AY k+l \ 2 ) - |E, fc (AY fc+1 )| 2 ]. 
1 Jt k 
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We can write E tfc |Y s - Y» +i | 2 < 2E tfc \Y S - Y k+1 \ 2 + 2E tk | AY k+1 \ 2 . By doing the 
same for X s — i , and taking 7 = C, we obtain 

|AF fc | 2 <(1 + Ch)E tk \AY k+1 \ 2 + Ch\AX k \ 2 + ChE tk \AY k+1 \ 2 

+ C[h 2 + E tk [ tk+1 \X S -X tk \ 2 ds + E tk [ tk+1 \Y s -Y tk+1 \ 2 ds] 
Jt k Jt k 

+ C[E tk [ tk+1 \Z s -Z tk \ 2 ds]. 

Using (9) yields |AF fc | 2 < (1 + Ch)E tk \ AY k+1 \ 2 + Ch\AX k \ 2 + Ch 2 . □ 
3.2 Proof o/E(E£o 1 ft k k+1 \Z» - Z t \ 2 dt)l = O(W). 
First of all, we can split this summation into two terms 

N-l ft ,,, N-l rt N-l 

E ( £ / K~Z t \Ht) P < CE( W \Z th -Z t \ 2 dt) P +CE(h £ |AZ fc | 2 ) P . 

A;=0 tk k=0 Jtk k=0 

Thanks to (9), we have E(j2 k =o l! k k+1 \Z tk -Z t \ 2 dtf < T p ~ l E k =o ti k +1 E|^ tfc - 
Z 4 | 2 frft = 0{hP). 

Scheme of the proof ofEfhY^ |AZ fc | 2 ) P = 0(h p ). The first key point is to 
slice the summation into small intervals and show that the result is true for 
small time intervals. The second key point is to use Rosenthal's inequality, see 
Theorem 2.12 page 23 of Hall and Heyde (if. By using (14) and taking the 
expectation, we can write : 

E(h ]T I AZ k \ 2 ) P <CE(J2 Var ifc AF fc+1 ) P + ChP. (15) 

k=0 k=0 

We use Rosenthal's inequality to upper bound 

fc1 fcl r \2p 

fe+1 7 



E( £ Var tfc AY k+1 ) P < CE( £ AY k+1 - E tk AY, 

k=0 k=0 

ki 

< C?, 2 ^ 1 \EAY k 2p +1 + EAY 2p + E(Y,{AY k - E tk AY k+1 



k=0 

By plugging this inequality into (15) and using the previous estimate on |AY/.|, 
we get 

E(h £ \AZ k \ 2 ) P < 0{hP) + CE(Y,{AY k - E tk AY k+1 )f P . (16) 

k=0 k=0 

We now tackle the term AY k - E tk AY k+1 . Using (11), we have £jj;Lo(A>A - 
^t k AY k+1 = Eto It k h+1 (^t k (f t N k - m)))ds. By doing the same kind of proof 
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as before, that is using the fact that / is Lipschitz and the results on E| AX k \ 2p 
and E\AY k \ 2 P, we find 

E( J^{AY k - E tk AY k+1 )f P < 0{hP) + C{hh) p E(h £ \AZ k \ 2 ) P . 

k=0 k=0 

By plugging this term back into (16), we can write (l— C (hki) p ^E(hYlk=o \AZ k \ 2 Y = 
0(h p ). Consequently, if we choose k\ < j— we come up 

withE(/i£fcLo|AZfc| 2 ) = 

(2C)Ph v ' 

0(h p ). This result can be extended to any summation involving at most Ak 
terms, where Ak < — ^— . We can cover the interval {0, • • • , iV — 1} with a fi- 

(2C)Ph 

nite number of elementary intervals of size Ak and we get E(ji J2^=o l^-^fcl 2 ) = 
0(h p ), which completes our proof. □ 
From this result and (9), we also deduce 

N-l 

E(h £ \AZ k \ 2 ) P = 0(h p ), (17) 

k=0 

which is very useful in the following. 



4 Proof of Theorem 7. 



To expand the error, we use usual techniques of stochastic analysis, combining 
martingale estimates and Malliavin calculus tools. 



4-1 Preliminary estimates 



Sections 4 and 5 contain proofs with similar calculations, which are quite 
technical. In order to be as clear as possible, we state two results really useful 
in the sequel, which are related to Malliavin calculus (see Nualart [l6(). The 
results give sufficient conditions for expectations and conditional expectations 
to be small w.r.t. the time step h. They are based on ideas from Kohatsu-Higa 
and Pettersson lB and Gobet and Munos 



Proposition 11 Let F G D 1 ' 2 with EijF| 2 +sup tfc < s < T E tk \D s F\ 2 < oo and let 
U be an It process of the form U t = U +fQa s ds+fQP s dW s , withsup tk<s<T E tk \a s \ 2 + 
su Pt.<s<T^tJAs| 2 < 00 ■ Then, V(t, t') such that t k < t < t' < t k+1 , 



\E tk [F(U t -U t ,)]\<(t'-t) (E tk \F 



]2( sup E tk \a s 

t<s<t' 

1 



+ ( sup E tk \D s F\ A )2( sup Etjfl 



t<s<t' 



t<s<t' 
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This proposition can be easily proved. Assume without loss of generality that 
F and U are one- dimensional. From the duality formula, we have E, tk [F(f^ a s ds+ 
ll (3 s dW s )\ = Et fc [// (Fa s +D S F • P s )ds\. Thanks to Cauchy Schwarz inequality 
and hypotheses on a and f3, we get the result. 

Definition 12 F satisfies the condition R k if F G and if C kyP (F) : = 

\\F\\ Lp +Ej<fcSup < ai) ... i8 .< T \\D 8u „ t8 .F\\ Lp < oo. 

Proposition 13 Let F satisfy the condition R 2 - For simplicity we set dW® = 
ds. Assume that U t G M d satisfies the following stochastic expansion property 

u t =t <S°(*) f , c5 2 (r)^;)^i, (P) 

i,j=0 Jo Jr >W 

where {(c^j 1 (t)) t > : < i, j < q, < i\ < 2} are adapted processes satisfying 

• V(i,j'),l < i,j < g, Vt G [0,T] ; c^°(t) satisfies R 2 , and 
C% P ■= sup < t < T sup 1 < iii < g C 2 , p (c^ (t)) < oo ; p > 1. 

• V(i,j), 1 < i, j < 9, Vt G [0,T] ; ^(t),^' ^),^^),^ 1 ^) safe/?/ and 
CYj, ■= sup < t < T supi^jx JCi iP (c5 1 (t))+Ci, p (c^(t))+Ci iP (c^°(t))+Ci iP (c^ 1 (t))} 
< oo ; p > 1. 

• V(i,j),0<i,j <q, Vt G [0,T], c§ 2 (t), cg(t), <§?(*) satisfy R , and C" p : = 
sup <,< T sup < M .< 9 {C , p (c^ 2 (t)) +C , p (c i ;' 1 (t)) +C 0lP (cgb°(*))} < oo ; p > 1. 

TTras, tnere zs a constant K{T) which depends polynomially on C2, P (F), C 2p , C^ p , C^ p 
(/or somep> 1) such that \E[FU t ] \ < K(T)h. 

Indeed, we have 

E(FU t ) = t E(FcS°(0 fc u i: }{s){f S c u i: ?{r)dWi)dWl) 



i,3=0 
1 r t rs 



Then, the result readily follows. 



i,j=l JO Jt i( s ) 

+ ttr E [ D i{Fcoj{t)}co']{s)co;f(r) 

j=l J Jr)(s) 

±f f]^[Dl{Fc^{t)c^{s)}c^{r) 

i=1 JO Jr)(s) 

+ f r EKwi^icjbV))* ds. 



cff(r) )dr ds 



dr ds 
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Remark 14 Under Hypothesis 2, we can show (see later the proof of (41) ) 
that for each t, X t — X t satisfies the expansion V. Hence, if F satisfies R2, 
Proposition 13 yields 

\nF(X t N -X t )]\ = 0(h) 
uniformly in t e [0, T], which is a very useful result for the sequel. 

4-2 Expansion ofY t ^ — Y tk 

In the following, we assume that Hypothesis 2 is in force. This implies in 
particular that u is bounded, of class C 3 ^ 2 ' 3 (see Theorem 4). We also easily 
prove that Vp > 1, V7c 6 {0, ■ ■ • , iV — 1} (see Nualart Q e.g.) 

• E tfc (sup \X t \ 2p ) <K(T)(l + \X tk \ 2p ), sup E tk ( sup \D s X t \ p ) < C, 

tk<t<T t k <s<T t k <t<T 

sup Et k ( sup \D r D s X t \P)+ sup E tk ( sup \D v D r D s X t \ p ) < C, 

t k <s,r<T t k <t<T t k <s,r,v<T t k <t<T 

(18) 

• Et k { sup \X t N \ 2p ) < K(T)(1 + \X"\ 2p ), sup E tk ( sup \D S X?\ P ) < C, 

t k <t<T N,t k <s<T t k <t<T 

sup Et fc ( sup \D r D 8 Xf \ p ) + sup E th ( sup \D v D r D s X t N \ p ) < C. 

N,t k <s,r<T t k <t<T N,t k <s,r,v<T t k <t<T 

(19) 

Due to the Markov property of (X^) k , one has Y/* = u N (t k , X^) for some 
Lipschitz function u N (t k ,-) (see Gobet et al. 0) with an obvious definition 
of u N . Actually, under our assumptions, this function is even three times 
differentiable w.r.t. x. Thus, the difference AY*, can be written as follows: 

AY k = (u N (t k ,X t N h )-u(t k ,X t N k )) + (u(t k ,X t N k )-u(t k ,X tk )). 

Since u is of class C^ 2 ' 3 , the last term of the previous inequality becomes 

u(t h ,X t N k )-u(t k ,X tk ) = V x u(tk,X tk )AX k + 0(\AX k \ 2 ). (20) 

To complete the proof, we apply the following lemma 

Lemma 15 Under Hypothesis 2, \u N (t k ,x) —u(t k ,x)\ < K(T,x)h. 

The result above is new but not so surprising. Indeed, if / is identically zero, 
the difference is only related to the weak approximation of $(X^) by $(X^): 
from Bally and Talay 0], one knows that this is of order h. 

The rest of this section is devoted to the proof of the lemma. We only give the 
proof for t k — 0. We want to find a upper bound for \u (0, x)—u(0, x) \ = \AY \. 

For the sake of clarity, we split the proof into several steps. 
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Step 1 : linearization of the error. We show that 

AY k = E tk (AY k+1 £ k + hf x {6 tk )AX k + h Xk ), (21) 

with 

Z k = (l + hf>(6 tk ) + f z (e tk )AW k ), (22) 



'y 

f tk+ \G (s,X s ) + f y {e tk )G y {s,X s ) + f' z (B tk )G z (s,X s ))ds (23) 

+ f(l - A) [AXidOtiAX,, + CK)( Y tk +1 ~ Y tk? + AZ k f';M k )&Z* k 
+ 2AX* k fl y {el){Y t N k+i - Y tk ) + 2AX*£«)AZ* + 2(Y t f +i - F tfc )4«)AZ* 

where #4 = X{t k ,X^,Y^ +x ,Z^) + (1 - A)0 tfc and G 0} G y ,G z are bounded 
functions. From (11) and by introducing f{0 tk ), we have 

AF fc = E tk (AY k+1 + h(f t N k - f(6 tk )) + J t [ k+ \f(0tk) ~ m))ds). (24) 

By applying It's formula to f(6 u ) between t k and s we show that, under 
Hypothesis 2 , E tk (f(6 tk )-f(6 s ))ds = h £? tjk (G (s,X a ))ds, where G 
is a bounded function. In the second term, perform a second order expansion 
of / around 9 tk to get 



Xk 



dX, 



fg - f(K) = fLVtJ&Xk + m k )AY k+1 + f'M k )AZ* k + fl(e tk )(Y tk+1 - Y t 



f k / 



+ J Q \l - A) [AX* k fL(0i)AX k + />, A J0C - Y tk f + AZ k f':M h )AZ* k 

(25) 



+ ZAXtf xy {9$ k ){Y» +i - Y tk ) + 2AXlf x M k )AZ* k + 2(Y» +i - Y tk )f y M k )AZ* k 

Note that E tk (Y tk+1 - F t J = E tk G y (s,X s )ds. If we closely look at (25), 
we can see that we need to develop AZ k . By using (5), we can write 



dX. 



Z t N k = -E tk (AY k+1 AW* k ) + -E tk (u(t k+1 ,X tk+1 )AWZ). 



±E tk (AY k+1 AW k *)+± 

Introducing the weak derivative of X tk+1 (see Nualart p. 109), the second 
term of this summation equals j^E tk J tk +1 V x u(t k+ i, X tk+1 )D t X tk+1 dt, where 
D t X tk+1 = V x X tk+r (V x X t )~ l a(t,X t ). Since Z tk = V x u(t k , X tk )cr(t k , X tk ), one 
gets 

AZ k = ^-E tk {AY k+1 AW* k ) 

+ - Jf E tfc (V^Ctfc+i, X tk+1 ) V x X tk+1 {V x X t )- l a{t, X t ) - V x u(t k , X tk )a(t k , X tk ))dt. 

The term in the second conditional expectation is equal to V x u(t k+ i, X t ) 
V x X tk+1 (V x X t )~ V(t, X t )±V x u(t, X t )a(t, X t )-V x u{t k , X tk )a(t k , X tk ): hence, 
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two applications of It's formula (for the first contribution between t and t k +i, 
for the second one between t k and t) prove that 



AZ* k = J* k+1 E tk (G z (s,X s ))ds + ^E tk (AY k+1 AW k ), (26) 

for a bounded function G z . Plugging this equality and (25) into (24) yields 
(21). 

Step 2 : another formula of AYq. First of all, we replace Y t ^ — Y tk by 
AY k+ i+Y tk+1 — Y tk in the expression of Xk- Then, easy computations combining 
Proposition 10 and estimates (9) show that 

Xk = IMS*) = O k (h) + 0(|AX fc | 2 + \AZ k \ 2 ). (27) 

From (21), we deduce the following equality 

AY = E(AY N £ • • • Cn-i + hY^(f' x (e ti )&Xi + X*)£o • ■ ■ (28) 

i=0 

Now it is enough to show that all terms of this summation are 0(h). In the 
following, i] = 1 and rji — £o • • • &-i f° r i < N. 

Step 3 : some results on t]n = £o • • • £jv-i> 
We establish the following results on r]N' 

r\ k satisfies the condition R 2 uniformly in k, i.e. Vfc, % e © 2 '°° 
and maxC 2 r>(%) < oo,Vp > 1, (29) 

k<N 

E( max |r7fe| p ) + supE( max |-D r r7fc| p ) + sup E( max \D r D s ri k \ p ) < oo. (30) 

Proof of (29). We have r/ = 1, and for i > 1 

77, = 77,^(1 + + (31) 

We begin to show that max fe <jv ||%||l p = 0(1) for p > 1. Since and 
are bounded, we easily prove that E ti-1 (l + ^(^.J + /^.JAWi-i)^ < 
(1+C/i), whence E|r/;| 2p < (l+C7i)E|T/;_i| 2p . We deduce that max fc <jv ||%||l = 
0(1). 

Now, let us show that max k <N E|-D r r/fe| p — 0(1), uniformly in r. Let r be such 
that tk-i < r < t k . Vi < fc — l : D r rji = 0. We note that -D,^ = rj k _if' z (6 tkl ). 
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For i > k + 1, we have 

D rVl = D rVi ^ + hD^J^)) + ^(rii-i^i^Afti, 

1=1 

- + ^E^fe/;(^)) +ee^(%/;^)) a ^- ( 32 ) 

j=k 1=1 j=k 

Applying Burkholder-Davis-Gundy's inequality to the martingale 
EUA-^C^AH/j yields 

E|D r ^| p < CTEWfc-i \ p + C p h £ El^^/i^.)) \ p + C ]T E|/i £ 

< CE|%_i| p + C/i^E|Z) r (^^(^))| p + C^2hY^E\D r (T] j f' Zi (6 t3 ))\ p 

j=k 1=1 j=k 

i-1 

KCil + Elrik-^+Ch £ E l^r> 

using the boundedness of the derivatives of /, maxj<jv = 0(1), idendity 
(7), u, cr e C^' 2 , and estimates (18). By applying Gronwall's lemma, we get 
maxfe^jvElA.^ < C(l + E|?7 fc _i| p ), t k _ x <r <t k . 

Then, maxfc<jv E\D r 7]k\ p = 0(1), uniformly in r G [0,T]. The proof concerning 
the derivative of order 2 can be done following the same scheme. □ 

Proof of (30). We begin to show that E(max fc <7v \i]k\ p ) < °o- The idea is 
to use a martingale property in order to apply Doob's inequality. Since rji = 

^-i/^^.JAW^l Thus, 



E(max < C(l + E(£ hfa-iWfytf^W + E(max | £ t^/^JAW^ 

— i=l — i=l 

The last term is upper bounded by CE(hT,f =1 h-ifi^u^W) $ < Ch 
J2iLi E|^i-i/^(^ i _ 1 )| p - Using the estimate (29), we get E(maxfc<7v \Vk\ p ) < oo. 
To prove that sup r<T E(max fe <7v \D r r] k \ p ) < oo, we proceed in the same way, 
by starting from (32). For the second derivative, this is analogous. 

Step 4 : we prove that E(AY N r) N ) = 0(h). 

If t}n were equal to 1, the results of Bally and Talay [l| would directly apply. 
Here the approach has to be different and we use techniques of Malliavin 
calculus. We have E(AY N rj N ) = E(r] N ^(X^) — r] N &(X T )). Let us introduce 
X t N ' x = (1 - X)X t + XX t N . Thus, we have 

E(AY nVn ) = f E(ri N & x {X%' x ){XZ - X T ))d\. 
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As $ G C 3+a , by using (29), (18) and (19), we note that r) N ^' x (X^ x ) satisfies 
R 2 - By applying Remark 14, we deduce that E(AY N r] N ) = 0(h). 

Step 5 : we prove that E(/^(# t JAXj?7j) = 0(h). This is a very similar proof 
to Step 4, in a case where <&(x) = x. 

Conclusion. We now work on hE(Y / ^ j 1 XiVi)-, where \xk\ < h+K(T, x)|AX fc | 2 + 
K(T,x)\AZ k \ 2 . Hence, 

N-l N-l N-l 

\h HXiVi)\ <C E H^h\)h 2 + K(T,x) £ ^(l^KIAX.I 2 + |AZ,| 2 )) 

i=0 i=0 i=0 

N-l 

<K(T,x)h + K(T,x) £ /iE(|^||A^| 2 ) 

i=0 

i iV-l i 

<K(T, x)h + K(T, x) (E( max t |i*|) 2 ) 5 (E(/i £ I a ^| 2 ) 2 ) 3 • 

— j=0 

By using (30) on (77^)^ and the upper bound (17) we get that \hE(Y / ^ l 1 XiVi)\ — 
K(T,x)h. By combining this result and the results of Step 4 and Step 5, (28) 
shows that |AY | < K(T,x)h. Lemma 15 is proved. □ 



5 Proof of Theorem 8. 



As it could be expected, its proof is more difficult. The main extra ingredient 
is the convergence of the weak derivative of the discrete BSDE (Y N , Z N ), with 
the rate of convergence N^ 1 ^ 2 . The next paragraph is aimed at proving this 
result. In the following, Hypothesis 3 is in force. 



5. 1 Proof of an intermediate result 

Proposition 16 Letr e]0,£i[. Under Hypothesis 3, we have max 1 ^ i ^^K\D r AY i \ 2 + 
/^(E^I 1 \D r AZ*\ 2 ) = 0(h), uniformly in r. 

This proposition is analogous to Theorem 6, where q — 2, and the scheme of 
its proof as well. However, there is a significative difference: the BSDE solved 
by the weak derivatives (see (33-34-35)) has a non Lipschitz driver, which 
requires extra technicalities that we detail. In what follows, we fix r g]0, ii[ 
and introduce some specific notations. X t stands for D r X t . In the case of Z tl 
which is a row vector, Z t is a matrix whose the i-th column is D l r Z^. It is 
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well-known (Proposition 5.3 of El Karoui et al. 0) that (Y t , Z t ) r <t<T solves 

% = & X (X T )T T + j\f x (9 s )X s + f y (6 s )Y s + f z (6 s )Z s )ds - (J t T z s *dw s y. 

(33) 

Regarding (Y N , Z N ), one obtains 

y? k =^[*K + ™*fZK + hV y f^ + hV t f»Zg\, (34) 
iJ^EUAW^j, (35) 

where we set V x f£ = V ' x j '(t k , X£,Y£ +i , Z£) and analogously for V y /£ and 
V z fu- Indeed, we can start from (4-5) and interchange conditional expecta- 
tions and weak derivatives (see Proposition 1.2.4 in Nualart (l6[). Another 
way to get (34-35) is to take advantage of the Markov structure of (X^) k to 
write Y t = y N (t k ,X^), where the function y is the solution of a dynamic 
programming equation, and then apply the chain rule. We omit further details. 

From (7), we also have 

Y t = V x u(t,X t )X t , % = V x (V x ua)*(t,X t )X t . (36) 

For the sake of clarity, let us write, for any process V, AV k = D r V^ — D r V th . 

't k , 



In particular, we have AZ k = D r (ZF* — Z* tk ) = Zj^ — Z tk , where Z tk is defined 



as hZ tk = E tfe Jt k +1 Z s ds (see the beginning of Section 3). 



5.1.1 Preparatory estimates 

In this part we give some L p -estimates (p > 1), which are repeatedly used 
inthe following calculations. 

. sup (E u \X»\ 2p )<C\Xp\ 2p , (37) 

i<3<N 

. E( ma <X jA^n = 0(1), (38) 

. V j e 0..N - 1, \Y~» | 2 < 0|xj| 2 , E( max \Yg \ 2p ) = 0(1), (39) 

• E(sup \X t \ 2p + sup |l^| 2p + sup \%\ 2p ) = 0(1), (40) 

0<t<T 0<t<T 0<t<T 

• Let F satisfy R 3 . Then, \R(F(X? - X t ))\ = 0(h). Furthermore, 

sup E|AA^| 2p = 0(h p ). (41) 

0<k<N 

• Analogously to (9), Vs G [tk,tk+i], we have 

E tk (\X S - JQ k \ 2p + \Y S - Y k \ 2p + \Z S -% k \ 2p ) = O k (h p ). (42) 
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Note that Xg = <r(0, x), and X t N k+i = (l+hb' x (t k: X^^U^jUh, X%)AWl)X t N k 

for 1 < k < N. Thus, we easily get E u \X^\ 2p < (1 + Ch)E u \X^ 1 | 2p , and (37) 
follows. The proof of (38) can be done as the proof of (30). 

Proof of (39). From (34), we use Young's inequality and boundedness of V/ 
to get 

j 2 < (1 + t^IE^^I 2 + C/.C/i + ^)(|^| 2 + ^Jl^J 2 + |^| 2 ). (43) 

From (35) and the Cauchy Schwarz inequality, we obtain h\Z£ | 2 < C(E ti \Y t f +1 \ 2 — 
|E t .F t ^ 1 1 2 ). Hence, with an appropriate choice of 7, (43) is reduced to \Y t f\ 2 < 
(1 + Ch)E t .\Y£ +1 \ 2 + Ch\X%\ 2 , and thus Gronwall's lemma yields 

JV-l 

\Y»\ 2 <CE u (\Y t N N \ 2 + h £ \X^)<C sup E ti \X»\ 2 . 

j=i i<j<N-l 

Finally, estimates (37) and (38) complete the proof. 

Proof of (40). E(sup < t < T \X t \ 2p ) = 0(1) follows from (18). The other esti- 
mates come from this result and (36). 

Proof of (41). Let us introduce X[ = V x X t (V x X r )~ 1 a(0, x) and write X^ — 
X t = X? -X' t + X' t -X t . 

Since X t = V x X t (V x X r )~ 1 a(r, X^, a direct application of Proposition 11 
with U t = cr(t,X t ) gives E(F(X' t -X t )) = 0(h) for F satisfying R 2 . Moreover, 
simple increment estimates yield sup i<T E|X t ' — X t \ 2p = 0{h p ). 

It remains to study the impact of the difference X™ — X' t . (X^) t > r and (X' t ) t > r 
are solutions of 

X?=a(0,x) + fVM*),X%s))X&)d8 + £, [\<T i )'Ms),x^ s) )x^ ) dw:, 

Jr i=l Jr 

X[ =(7(0, x) + fb' x {s,X s )X l s ds + Y. f\°i)' x {s,Xs)X' s dWl. (44) 

Jr i=1 Jr 

For the sake of simplicity, we take 6 = and d — q — 1. If we set cr'(s) = 
Jo a 'x( s i X s + A(Xj v — X s ))dX, we observe that AX t solves the linear equation 
AX t = ft[<j(r](s),X% s) ) - a(s, X?)]dW 8 + Jo a'(s)AX s dW s , which solution is 
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given by (see Theorem 56 p. 271 in Protter |lj 

AX t =e t j\^[a(v(s),Xl ) ) - a(s, X»))(dW s - a'(s)ds) 

+ (a' t (v,X») + ^ x (v,X»)a\ V (v),X» (v) ))dv](dW s - a'(s)ds) 
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where e t = 1 + Jq a'(s)e s dW s . This proves that AX t satisfies the property V. 
Analogously, if we define c"(s) = f a xx {s,X s + A(Xj v — X s ))d\ and ef = 
1 + f* <j' x (s, Xf)e^ dW s , simple computations lead to 

X? - A =<? fa)-\teWs),X%.))>&) - o> x (8,X?)X»] + a"(s)X' s AX s ) 
\dW s -a x (s,X?)ds). 

From the above representation, it is straightforward to conclude sup t<T E| AX t \ 2p = 

0{h p ). Now, let us upper bound K(F(X^ — X' t )) which can be decomposed 
into several terms. 

. The contribution associated to ef £(ef ^XM*), ^J*))^-^ X»)X?](dW s 
a' x (s, X^)ds) satisfies property V, thus Proposition 13 yields the expected 
result. 

• The contribution E(Fef J r *(ef )- 1 a"(s)X' s AX s a' x (s, X? )ds) is equal to 
X? E(Fef (ef )- 1 a"(s)X' s AX s a' x (s, xf))ds = 0(h) in view of Remark 14. 

• In the same way, the duality relationship ensures that the last contribution 
E(Fef J*(e^a"(s)X' s AX s dW s ) = f* E(D s (Fe^)(ef)~ 1 a"(s)X' s AX s )ds is 
& 0(h) (using here that F satisfies R 3 ). 

Proof of (4^)- In view of X t = D r X t = x X t (^7 x X r )~ l a(r, X r ), the estimate 
on the increments of X t becomes clear. The other ones easily follow. □ 

5.1.2 Proof o/maxi^^ElAFil 2 = 0(h). 

Assume that for some non negative random variable = O k (h) + |AXfc| 2 + 
\AZ k \ 2 , one has 

\AY k \ 2 < (l + Ch)E t JAY k ~^ 1 \ 2 + h\AX k \ 2 + hA k O k (l). (45) 

Take the expectation on both sides, use estimates (41) and those of Proposition 
5 to get 

N-l 

E\AY k \ 2 < CE|A1^| 2 + 0(h) + ChJ2 E(|AZ fc | 2 O fc (l)). 

fc=0 
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On the one hand, as AY/v = &(X")X t N N - &(X t „)X tin clearly E\AY^\ 2 = 
0(h). On the other hand, in view of (17) with p = 2, the summation above is 
a 0(h). This proves maxi< fc <Ar E| AY k \ 2 = 0(h). 

Proof of (45). From (33) and (34), we obtain 

A? fe = E tk (Ai£i) + ®t k ( £ +1 \y x f£X? h - f'M)X s 

Since / G C 2 ' 4 ' 4 ' 4 , it follows that for any 7 > (to be fixed later) 

|AF fc | 2 <(l+ 1 h)\E tk (AY k+1 )\ 2 + C(h+-)E tk ( / [\V x f t N k X t N k -f' x (9 s W 

+ \^vft N A - m)Y s \ 2 + |V,/^ - £(*.)5|>) (46) 
<(1 + 7/^(^0 1 2 + C(Zi + l -)(Tl + T 2 ), (47) 

where we put T k = E tk (J^[\x][ - X s \ 2 + |^ - Z\ 2 + K " 

^ = /i fc+i (^+i^-^iri 2 +i^— v;^, i 2 +i^ s — ^i 2 )(|j^;i 2 +i^i 2 -hi^;i 2 )^. 

To get (45), we need to simplify (47), by estimating T k and T 2 . 

Term T\. Firstly, wewrite E t J>^-^| 2 < 2Ej>^->;| 2 + 2E,JAlCi| 2 - 

We do the same for Xf k — X s . Then, the usual increment estimates yield 

Et fc - ^| 2 + E t jxJ - X s \ 2 < O k (h) + 2\AX k \ 2 + 2E tk |AnTi| 2 - 
Secondly, analogously to (13), we have 

E, fc / \Z t N k - Z s \ 2 ds = E tk / \Z tk - Z s \ 2 ds + hE tk \Z» - Z t f. 

Finally, we obtain T 1 < Ch(O k (h) + |AX^| 2 + E t JAY^i| 2 + |AZ fc | 2 ). 

Term T) 2 . Easy calculations combining (9), Proposition 10 and (40) give T k < 
(O k (h 2 ) + h\AX k \ 2 + h\AZ k \ 2 )O k (l) = hk k O k (l). 

Conclusion. Plugging the estimates on T k and T k into (47), we get 

|A? fc | 2 <(1 + 7 ^)|E, fc (AlCi)| 2 + Ch(h + -)|A^| 2 

7 

+ C^ + i)(|^|- + ^|^|- + AA(l))^ (48) 

Note thatjg = E tk (AW k (Y^ + f t k k +1 [f x (0 s )X s + f'(9 s )Y s + f' z (9 s )Z s ]ds)), 
whence hAZ tk = E tk (AW k (AY k ~^ 1 + [f^(9 s )X s + f y {6 s )Y s + f z {6 s )Z s ]ds)) . 
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By proceeding as before, we easily prove 

h\A^Z\ k \ 2 < C(E t JA>£| 2 - lE^Aldl 2 ) + O k (h 2 ). (49) 
Combining this upper bound with (48) for a good choice of 7 gives (45). □ 



5.1.3 Proof of hE ( J2 k =i I Az k 1 2 ) = O (h) . 

In view of (42), this is equivalent to prove hE^k^i |AZ fe | 2 ) = 0(h). To 
establish this estimate, we start from (49) to get 

N-l . N-l 

h £ E|AZ fc | 2 <CJ2 (E|A? fe | 2 - E|E tfc Ay fc+1 | 2 ) + CE|A1^| 2 + 0(h). (50) 
fc=i fc=i 



Now, we work on |AF fe | 2 - |Ei fc AF fc+ i| 2 . The choice 7 = 2C 2 in (48) leads to 

|An| 2 -|E tfc (An7i)| 2 < 7* fc (A>Ci)| 2 + h(^ + C/i)|A^| 2 
+ KCh + ^)(|a3g| 2 + E t JAn^| 2 + A fc O fc (l)). 

From (41) and the result from Section 5.1.2, we have maxi<jt<jv E(| AX^] 2 + 
|AF fc | 2 ) = 0(h). We also have E(A fc O fc (l)) = 0(h) + E(| AZ k \ 2 O k (l)). Conse- 
quently, for h small enough, one has E|AF fe | 2 -ElE^Al^)! 2 < ^E\A~Z k \ 2 + 
0(h 2 ) + ChE(\AZ k \ 2 O k (l)). Putting this estimate into (50) yields 

1 N-l _ N-l 

-h ]T mz k \ 2 < O(h) + C7i ]T E(|AZ fc | 2 O fc (l)). 

6 k=l k=l 

Inequality (17) with p = 2 directly shows that the sum above is a 0(h). □ 



5. 2 Expansion of Z\° — Z tk 



We recall that u G C 2 ' 4 owing to Hypothesis 3. From (26), we have AZ k = 
0(h) + ±E tk [(u N (t k+1: X t N k+i ) -u(t k+1 ,X tk+1 ))AW£]. Let (Xf) t > s denote the 
solution of the SDE (1) starting at time s from x. We write X t for X®' x . Note 

t Xt 

that X tk+1 = X tk ' +1 tk . In the same way, the Euler scheme starting at time t k 
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at x is denoted by (X^ tk,x )j> k . With this notation we can rewrite AZ k 

AZ fe = x E,J(^(t fc+1 ,X, fc ; i ' tfc ) -u(t k+1 ,X tk+1 ))AW*} + 0(h), 

1 t 

=-E, fc [( M (t fc+1) <; i tfc )- M (4 +l! X tfc+1 ))A^] 

1 Nt x N t x N 

+ -E tk [(u N (t k+1 ,X tk ' + k 1 ' tk )-u(t k+1 ,X t k k+i tk ))AW* k } + 0(h). (51) 

We work on the first two terms separately by proving 

t x N 

Lemma 17 \E tk [( u (t fc+1 , -«(**+!, X tfc+1 ))A^ fc *] = 0(| AX fc | 2 )+0(^) 

+ [V :c (V :cM a)*(t fe ,XiJAX fc ]*. 

TV £ X^ t X^ 

Lemma 18 1^ [K(4 +1 , X^' *) - x4*'=))AW*]| = O k (h). 

The combination of these Lemmas completes the proof of Theorem 8. 

5. 2. 1 Proof of Lemma 1 7. 

For the sake of simplicity, let A N X k+ i denote X tk+l k —X tk+1 (which is different 



-N,t k ,X" 



from AX k+ i = X tk+1 k —X tk+1 ). From a Taylor-Lagrange formula, we obtain 

t k x N 

u(t k+1 ,X tk+l k ) - u(t k+1 , X tk+1 ) = u x (t k+1 , X tk+1 )A N X k+1 
+ f {1 - X)(A N X k+1 )*H x (u)(t k+1 ,X tk+1 + XA N X k+1 )A N X k+1 d\. 

J 

Thus, using the duality relationship, one has 

E tk [(u(t k+1 ,X t k k+1 tk ) - u(t k+1 , X tk+1 ))AW* k ] 

= / R\(t)dt + / R 2 k (t)dt + / (1 - X)R 3 k (X)dX, 
Jt k Jt k Jo 

with Rl(t) = E tk [(A N X k+i yH x (u)(t k+1 ,X tk+1 )D t X tk+1 }, 
R 2 k (t) = E tk [u' x (t k+1 ,X tk+1 )D t (A N X k+1 )}, 

R 3 k (X) = E tk [(A N X k+1 )*H x (u)(t k+1 ,X tk+l + XA N X k+1 )A N X k+1 AW:\. 
Expansion of R k (t). Clearly A N X k+ i = AX k + U tk+1 —U tk , where U is an Ito 

tfc }X^ 

process with drift term a s = b(s,X s k ) — b(s,X s ) and diffusion term (3 S = 

tfc jX^ 

a(s, X s k ) — a(s, X s ), both being bounded. Thus, we can apply Proposition 
11, letting F = H x (u)(t k+ i, X tk )D t X tk . Because u G C^' 4 and in view of 
(18), we get 

Rl(t) = 0(h) + (AX k )*E tk [H x (u)(t k+1 ,X tk+l )D t X tk+l \. 
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We expand the latter factor. As D t X tk+1 = V x X tk+1 (V x X t ) X t ), we have 

H x (u)(t k+1 ,X tk+1 )D t X tk+1 = (H x (u)(t k+1: X tk+1 )<j(t,Xt) - H x (u)(t,X t )a(t,X t )) 

+ (H x (u)(t, X t )a(t, X t ) - H x (u)(t k , X tk )a(t k , X tk )) 
+ (H x {u){t k+1 ,X tk+1 )[V x X tk+1 {V x X t )-' - I]a(t,X t )) 
+ H x (u)(t k ,X tk )a(t k ,X tk ). 

The first three contributions in the r.h.s. above can be handled in the same 
way and we give a detailed proof only for the first one. It is enough to apply 
Proposition 11 with F = a(t,X t ) and U s = H x (u)(s,X s ). Then, E tk [F(U tk+1 - 
U t )} is of order h with a constant involving b,a,u and its derivatives up to 
order 4. Finally, this gives 

R\{t) = 0(h) + (AX k )*H x (u)(t k ,X tk )a(t k ,X tk ), 

uniformly in t G [t k , t k+1 }. 

Expansion of R\(t). For t k < t < t k+1 , we have 

vJV + vN + vN + 

D t (A N X k+1 ) = [V,X^* k (V x X t v Y 1 - I]a(t,X t ^ k ) 

- [V x X tk+1 (V x X t )' l -I]a(t,X t ) - (a(t,X t )-a(t k: X tk )) 
x N t 

+ a(t,X t k ) - a(t k ,X t N k ) + a(t k ,X») - a(t k ,X tk ). 

As before, apply Proposition 11 to each of these terms but the last one, with 
F = u' x (t k+ i, X tk+l ), using u,b,a G C^' 4 and (18). It follows that R k (t) = 
0(h) + E tk [u' x (t k+1 ,X tk+1 )](a{t k , X t N J - a(t k , X tk )). An application of It's for- 
mula yields 

Rl(t) = 0(h) + j2u' x St k ,X tk )(a\t k ,X? k )-a\t k ,X tk )) 
i=i 

= 0(h+\AX k \ 2 ) + iru' x V x ([ar)(t k ,X tk )AX k , 
i=i 

uniformly in t G [t k ,t k+ i\. Finally, simple matrix computations lead to 
Rl(t) + R 2 k (t) = 0(h + |AX fe | 2 ) + [W x (W x ua)*(t k , X tk )AX k ]*. 



Upper bound for Rl(X). To complete the proof of Lemma 17, note that it 
remains to justify that R k (\) = hO(h + \AX k \ 2 ) uniformly in A. The duality 
formula gives 

R\(X) =E tk [f tk+1 D t [(A N X k+1 )*H x (u)(t k+1 ,X tk+1 + XA N X k+1 )A N X k+1 ]dt. 
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The term in the integral equals Y^tj=i[^D t (A N X k+lii )A N X k+1J dl ijX .u(t k+1 , X tk+1 + 
XA N X k+1 ) + A N X k+1>i A N X k+ltj b t (dl i!X .u(t k+1 ,X tk+1 + AAjvXfc+i))]. Thanks 
to (18) and (19) and successive applications of Proposition 11, we finally prove 
our assertion. We omit further details. □ 



5. 2. 2 Proof of Lemma 1 8 

As for Lemma 15, we only do the proof for t k = 0, i.e. we have to show 
\E tk X^ x )-u{h, X»?)) AW£] | < K(T, x)h\ We have E[(u N (t u X?**)- 

u(t u X^ x ))AW^] = E[AYiAWq]. By using (21), we come up with 

EIAFxAWJ] = E[£ 1 ..4 N - 1 AY N AW5] + E[h ^(f'M^AX, + ^i-^-iAW*}, 

i=i 

where Xi — ^uiXi) (& and Xi are defined in (22) and (23)). In the following fji 
denotes an d Vi = 1- We easily prove that (fji)i<i<N has the analogous 

properties to (r)i)o<i<N- Estimates (29) and (30) remain valid for fj and under 
Hypothesis 3, the estimate (29) becomes 

rf k satisfies R 3 uniformly in k. (52) 



Step 1 : Proof of E[^ 1 ...^ N _ 1 AY N AW*} = E[rf N AY N AW*} = 0(h 2 ). 
As before, we use the duality formula: 



E[7f N AY N AW*} =E f l (D t [rf N \AY N + rf N D t [AY N ])dt. 
Jo 

Since rfw satisfies (52), we proceed as in Step 4 of Lemma 15 and we get 
E(D t [r]~N] AYn) = 0(h). Furthermore, we have 

D t [AY N \ = - <$>'(X T ))D t X$ + $'(I T )(AIr " D t X T ). 

On the one hand, analogously to previous computations, we establish E(rf N (^'(X^)- 
&{X T ))D t X*!)=0{h). 

On the other hand, we prove E(rf N ® (X T )(D t X$ - D t X T )) = 0(h). Thanks 
to (18) and (29), rf N ^'(X T ) satisfies condition R 3 . Then, by applying (41), we 
get the result. 

Step 2 : Proof of E[/i£^V*(0*JA^ 1 ...&_ 1 AW*] = 0(h 2 ). 
This is a similar proof to the one done at Step 1, with = x. 

Step 3 : Proof of E[h EjlT 1 M*AW *] = 0(h 2 ). 

A careful inspection of the definition of Go, G y and G z appearing in (23) 
shows that under Hypothesis 3, these functions are continuously differentiable 
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w.r.t. the variable x (with a bounded derivative). Hence, if we write Xi = 
x! + /o(l-A)x?(A)dAwith(see(23)) 



x\ = 



[ U+1 (G (s,X s ) + f y (6 u )G y (s,X s ) + f' z (9 ti )G z (s,X s ))ds, 

ax* fLie^AX; + f"M)(Y". - k.) 2 + az f'!M.)Az 



xl(\) = axiC^ax, + f';M)( Y u + i ~ Y ^) 2 + azj'Uo^z* 

+ 2AX*/> t A J(y t f +1 - Y ti ) + 2AX*f: z (0*)AZ* + 2(F t f +i - ^(^JAZ*, 

we note that the random variable Xi is m IO 1 ' 00 - Thus and because x% — E^(x«), 
one has E[M*AW *] = E[ X ^AW*] = E^ 1 ( Xi + WtXi)^]. 

The upper bound X i = ®u(Xi) = O l (h) + 0(|AX,| 2 + |AZ,| 2 ) (see (27)) is 
sufficient to show EfX^ 1 XiDtVi] — 0(1) uniformly in t (follow the arguments 
of the conclusion of the proof of Lemma 15 and use (30) with rj). 

Now, it remains to establish EfE^ 1 JfcAXi] = 0(1). On the one hand, clearly 
E ti [AXi] = Oi(h) and we conclude Ef^^i 1 Vi^tXi] = 0(1) uniformly in t. On 
the other hand, xf can be decomposed into several contributions, which can be 
analyzed with the same arguments. Let us detail how to handle one of them, 
for instance EfE^ 1 fjiD t (AX* f'^JAZ*)] which has to be a 0(1). We do 
the proof for d = q = 1. Write D^AX^MJAZ,) = AXj'^MjD^AZ,) + 
A(AX i )/; z (^)AZ i + AX i A(/«(^))AZ i . As /" is bounded, we have 



JV-1 



JV-1 



|e[£ vAx^MJd^az^ < e[ £ I^I|axh/; 2 (^)||a(a^; 
i=i i=i 

JV-1 i JV-1 

< c(e( £ (l^n AX,| 2 ))) 5 (e( £ | A(AZ< 



2\\ 2 



1=1 



1=1 



Thanks to Proposition 16, (52) and Proposition 5, we get that 
EfE^T^AX/^XAA^)] = 0(1). Analogously, using (52-17-41), we ob- 
tain ElxZ^l 1 rjiiDtAX^f'^e^AZi] = 0(1). It remains to demonstrate that 

mil 1 vAXiDt^MMk = 0(1). We have 



Dt(LMJ) = f xz MX XD t x u + (i - 

+ C(^)(AA^ t f +1 + (1 - A)A^) + fZM^vK + (1 " A)A^). 



The most difficult term to bound among these three ones is the one which 
contains D t Zf_ . If we write \D t Zf. + (1 - \)D t Z t . = \D t (AZ t ) + D t Z u , we 
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obtain 

i=i 

N-l i N-l i 

<0(e(]T |A(A^)| 2 )) 5 (E(5:(|AX l | 2 |^| 2 |AZ J | 2 ))) 5 , 

i=l i=l 

N-l i N-l i i 

< c(E(£ |A(AZ i )| 2 )) 5 (E(E |A^| 2 ) 2 ) 3 (E(maK |^| 4 max |AX,| 4 )) 3 , 

i=l i=l ~ ~ ~ ~ 

Applying Proposition 16, (17), Proposition 5 and (30) (with fj) lead to 
E[E^I 1 »7iAX i /;;(^)AA(AZ i )AZ i ] = 0(1). Proposition 5, (17), (38), (39), 
(40) and (52) enable us to prove that the others terms of 
EfEill 1 *AX i A(^(fiJ))AZ i ] are 0(1). □ 
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